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ABSTRACT 

Recovering the Cosmic Microwave Background (CMB) from WMAP data requires galactic foreground emissions to be 
accurately separated out. Most component separation techniques rely on second order statistics such as Internal Linear 
Combination (ILC) techniques. In this paper, we present a new WMAP 9-year CMB map, with 15 arcmin resolution, 
which is reconstructed using a recently introduced sparse component separation technique, coined Local Generalized 
Morphological Component Analysis (LGMCA). LGMCA emphasizes on the sparsity of the components to be retrieved 
in the wavelet domain. We show that although derived from a radically different separation criterion (i.e. sparsity), the 
LGMCA- WMAP 9 map and its power spectrum are fully consistent with their more recent estimates from WMAP 9. 
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1. Introduction 

The Cosmic Microwave Background (CMB) is a snapshot 
of the state of the Universe at the time of recombination. It 
provides information about the primordial Universe and its 
evolution to the current state. Our current understanding of 
our Universe is heavily based on measurements of the CMB 
radiation. The statistical properties of CMB fluctuations 
depend on the primordial perturbations from which they 
arose, as well as on the subsequent evolution of the Universe 
as a whole. For cosmological models in which initial pertur- 
bations are of a Gaussian nature, the information carried by 
CMB anisotropics can completely be characterized by their 
angular power spectrum which depends on the cosmologi- 
cal parameters. This makes the precise measurement of the 
CMB power spectrum a gold mine for understanding and 
describing the Universe throughout its history. 

In the estimation of the CMB map, the astrophysical 
foreground emissions from our galaxy and the extragalactic 
sources have to be removed. In addition, the instrumental 
noise hinders the estimation of the CMB map. In the low 
frequency regime (below 100GHz, i.e. for WMAP chan- 
nels) the strongest contamination comes from the galactic 
synchrotron and free-free emission, with the highest contri- 
bution at large angular scales. At higher frequencies, dust 
emissions dominate whereas the synchrotron and free-free 
emissions are low. The spinning dust is an extra emission 
which spatially correlates with dust and dominates at low 
frequencies. 

Since second-order statistics provide sufficient statistics 
for a gaussian CMB field, most of the component separation 
techniques, such as the Internal Linear Combination (ILC), 
are built upon them to recover the CMB map from the ob- 
served sky maps. However, these techniques are not opti- 
mal for non-stationary and non- Gaussian components such 
as the foregrounds (or even non-stationary noise). On the 



contrary sparsity-based source separation techniques that 
focus on the higher-order sta tistics of the components have 
prove n to be highly efficient (Bobi n et al.|2 007; Bobin et al. 
20131). 
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In this paper, we present a new WMAP 9-year CMB 
estimation based on this sparsity concept, and we compare 
the results to the official WMAP products. Next section [2] 
briefly describes the LGMCA method. We then describe in 
section [3] the processing of WMAP data and the derived 
LGMCA products are displayed in section [4] Results on 
WMAP simulations are also presented in Appendix |C| 



2. Component Separation for CMB maps 

Exploiting the fact that foreground components are sparse 
in the wavelet domain (i.e. a few wavelet coefficients are 
enough t o represent most o f the energy of the component), 
LGMCA ( Bobin et al.|2 013) estimates both the components 
of interest and the mixing matrix by maximizing the level 
sparsity of each component; it seeks the sparsest sources 
possible in a wavelet basis. The assumption is that the ob- 
served sky is a linear combination of all the components, 
each resulting from a completely different physical pro- 
cess, and the instrumental noise. The separation principle 
in this method relies on the different spatial morphologies 
or structures of the various foregrounds, which translate 
into different sparsity patterns when transformed to a fixed 
wavelet dictionary. A linear combination of these compo- 
nents decreases the level of sparsity. Therefore, reconstruct- 
ing each component from the observed map, by maximizing 
its sparsity in wavelet space, is an efficient strategy to dis- 
tinguish between physically different sources. As mentioned 
in Appendix |Aj channels resolution variation and spatial 
variations of component emissions such as dust are taken 
into account by estimating a mixing matrix per wavelet 
scale and per area. Full details are given in |Bobin et al. 
([2013) 
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In Appendix [Cj LGMCA is evaluated on simulated 
WMAP data and it is found that: 

— The recovered power spectrum from the LGMCA map 
at 15 arcminutes is within the 2 a error bars from the 
input CMB power spectrum. 

— Propagated noise in the map is the main residual con- 
tamination for LGMCA ; for low multipoles, the con- 
taminants are much lower than cosmic variance. 

— Compared to a pixel-based localized ILC computed at 
1 degree, both noise and foregrounds residuals are lower 
using LGMCA. 

— No significant non-gaussianities at various scales and 
positions are detected in LGMCA maps, either at 1 
degree or 15 arcminutes. Compared to the pixel-based 
localized ILC, no significant difference is observed for 
LGMCA at 1 degree compared to the errors expected. 

Altogether, these results incentivize applying LGMCA to 
WMAP9 data. 



3. Map and Power Spectrum Estimation 



ILC estimate of the CMB map at 1 degree 



dust template is the most important one, as it helps in 
removing the dust emission on small scales, which is oth- 
erwise significant in the W channel. As the spinning dust 
is spatially correlated with thermal dust, the thermal dust 
template also helps in reducing the spinning dust residu- 
als. It is usually standard to u se the 408 MHz synchrotron 
map of (Haslam et al. 1981) as a template for the syn- 
chrotron emission. However, this map has quite a low reso- 
lution of about 1 degree. Furthermore, adding this template 
to LGMCA did not show any improvement for the compo- 
nent separation. 

LGMCA map: following QBobin et aTT] |2Q13[ ) ? the 
LGMCA is applied to the five offset-corrected WMAP maps 
and the two templates used as extra observations. Details 
on LGMCA parameters are given in Appendix [Bj 

LGMCA power spectrum: using the mixing matri- 
ces previously estimated from all data, we can estimate 
a CMB map from each of the nine individual year data 
set. From these nine maps, we have derived all 36 possible 
cross-spectra using the high resolution temperature anal- 
ysis mask kq85 with (f sky = 0.75)[] provided by WMAP 
collaboration. The final CMB power spectrum is then ob- 
tained by aver aging them and by using a MASTER mask 
deconvolution ( Hivon et al.||2QQ2 ). In contrast to the CMB 
and foreground signals, noise is uncorrelated between dif- 
ferent years of data and it is therefore considerably reduced 
in the averaged cross spectrum. 




Fig. 1. ILC official WMAP 9 years map (1 degree resolu- 
tion). Units in \iK. 

The WMAP satellite has observed the sky in five fre- 
quency bands denoted K, Ka, Q, V and W centered on 
23, 33, 41, 61 and 94GHz respectively. The released data 
includes sky maps obtained with ten differencing assem- 
blies, for nine individual years; per year, there is one map 
for the K band and one for the Ka band, two for the Q 
band, two for the V band and four for the W band. These 
maps are sampled using the HEALPix pixelization scheme 
at a resolution corresponding to nside of 1024. At first, 
we average all the differencing-assembly maps obtained for 
the same frequency band, which yields five band- averaged 
maps. However, these maps are not offset-corrected. In or- 
der to determine the offset value for a particular frequency 
band, we use the standard resolution of the 9-year band- 
average maps of WMAP as reference maps. Offset values 
are obtained by determining the mean of the difference be- 
tween the band-averaged map being considered, and the 
reference map. 

At the WMAP frequencies, the major sources of con- 
tamination in the maps are the synchrotron, free-free, spin- 
ning dust and thermal dust. To model the foreground 
contamination, we use two foreground template s in our 
analysis: du st at 100 microns, as obtained by (Schlegel 
et al. 1998[) and the composite all-sky H-alpha map of 
(|Fink 



4. Results 

LGMCA Map and Power Spectrum 

Figure [I] shows the official ILC WMAP 9 years CMB map, 
which has a 1 degree resolution, and Figure [2] shows the 
LGMCA WMAP CMB map at 15 arcmin resolution and 
at 1 degree. The 15 arcmin LGMCA map exhibits slight 
high-frequency structures which are likely related to point 
sources emission on the galactic center. The two 1 degree 
maps look very clean, even on the galactic center. Figure [3] 
features the difference at resolution 1 degree between the 
CMB map estimated by LGMCA and the official ILC-based 
CMB map provided by the WMAP consortium. The dif- 
ference map shows significant foreground residuals on the 
galactic center, that are however roughly three times lower 
than the CMB level. When masking the galactic center with 
the kq85 WMAP mask (Fsky=75%), no significant feature 
can be seen anymore. 

The estimated power spectrum is displayed in Figure 4] 
The errors bars, essentially the cosmic variance and the 
noise-related variance, are derived from c lassical power 
spect rum variance estimators as described in ( |Greason et al.| 
|2013[ ). Note that except for i < 32, the official WMAP 9 
years power spectrum (in blue in Fig. [4| is computed from 
the V and W bands only. In opposite, the spectrum in red 
in Fig. [4] is derived from the full dataset. Even if the esti- 
mation procedure differs slightly, it is remarkable that the 
two power spectra look very similar. They also tend to de- 
part from each other with a slightly lower third multipole. 
As well the measurement about the third peak seem to be 
higher to some extent. However both spectra are compati- 
ble at all scales within 2 a error bars. 

Sanity check: The LGMCA algorithm computes mixture 
parameters (i.e. mixing matrices and their inverse) which 



i>einer| 2003). Among these templates, the thermal 1 http:/ /lambda. gsfc.nasa.gov/product/map/dr5/m_products.cfm 
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LGMCA estimate of the CMB map at 15 arcmin 



LGMCA estimate of the CMB map at 1 degree 



Fig. 2. Estimated LGMCA CMB map from WMAP (9 years) at 15 arcmin resolution and 1 degree. Units in \xK. 



Difference between LGMCA and ILC map at 1 degree 



Difference between LGMCA and ILC map at 1 degree 




Fig. 3. Difference between the estimated CMB map with LGMCA and the official ILC map. The bottom panel shows 
the same difference map masked with the kq85 mask. 
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Fig. 4. Estimated CMB map power spectrum from WMAP (9 years) in linear scale (top panel) and logarithmic scale 
(bottom panel). Units in fiK 2 . 



Multipoles(l) 

Fig. 5. Mean ratio of Cj T /Cj T ' th over 100 random CMB 
simulations. The error bars are related to the cosmic vari- 
ance over the 100 simulations. 



are then applied to the data to estimate the CMB map. It is 
important to note that the final CMB map linearly depends 
on the input data. This makes it possible to check whether 
the inversion process may induce some bias at the level of 
the estimated CMB map and its power spectrum (assum- 
ing no calibration or beam errors). For that purpose, we 
applied exactly the same parameters we computed from the 
real data to 100 random CMB realizations. These realiza- 
tions were generated as Gaussian random processes with a 
power spectrum defined by the WMAP 9-year best-fit theo- 
retical power spectrum. The point spread function of these 
simulations are chosen as the beams of the 9-year band- 
average maps of WMAP thus providing pure CMB simu- 
lations that mimic the CMB part of the WMAP 9 years 
data. Figure [5] shows in blue the ratio between the esti- 
mated and theoretical power spectrum Cj T /Cj T,th com- 
puted from an average of 100 random CMB simulations. 
The theoretical power spectrum Cj T,th of the simulated 
CMB maps is exactly defined as the WMAP 9-year best- 
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Wavelet Scale Skewness 



■ ILC Skewness 
► LGMCA Skewness 



Wavelet Scale Kurtosis 



■ ILC Kurtosis 
► LGMCA Kurtosis 



Wavelet scale 



Fig. 6. Comparison of sk ewne ss and Kurtosis in ILC and LGMCA map at 1 degree computed for the various wavelet 
scales described in Figure C.4 (left). The wavelet filters peaks respectively about multipole 240, 120, 60, 30 and 15 for 
scales 1 to 5. 



fit theoretical power spectrum. The error bars are defined 
as the cosmic variance from 100 simulations. If the estima- 
tion by LGMCA of the CMB map and more specifically 
its power spectrum were biased, the Cj T /Cj T,th should 
depart from 1. Figure [5] shows that there is no statistical 
evidence of discrepancy from 1. We therefore conclude that 
the LGMCA does not introduce any bias at the level of the 
CMB power spectrum. Any bias of the estimated CMB 
map will come from remaining noise and foreground con- 
tamination. 

Higher Order Statistics — Non-Gaussianities: Higher order 
statistics were also computed for the full 9 year CMB maps 
recovered at 1 degree by ILC and LGMCA to assess poten- 
tial differences in their distribution. The 75% mask was em- 
ployed on both maps to avoid computing the higher order 
statistics in regions contaminated by foreground residuals. 
Sparse inpainting was then performed t o interpolate the 
signal inside the mask ( Starck et al.|2013 ) in order to avoid 
artifacts on the wavelet coefficients, and the skewness and 
kurtosis were calculated for each wavelet scale considering 
only wavelet coefficients within the mask. These statistics 
were centered and normalized by processing similarly a set 
of 100 realizations of noise and CMB according to the A- 
CDM fit provided by the WMAP consortium. Figures [6] 
shows the skewness and kurtosis versus the wavelet scale, 
and illustrates that not only both ILC and LGMCA maps 
are compatible with no non-gaussianities, but also that 
no significant difference between ILC and LGMCA can be 
found with these statistics at that resolution. 



5. Conclusion 

We have investigated how sparsity could be used for 
WMAP CMB map reconstruction. Based on WMAP sim- 
ulations, we have shown that LGMCA provides a low fore- 
ground map, and that noise remains the major source of 
contamination. Then a high resolution (15 arcmin) clean 
CMB map has been computed from the full WMAP 9 
years dataset, and its power spectrum was estimated. 
Remarkably, though LGMCA-based and official WMAP 9 
years power spectrum were derived from completely dif- 
ferent estimation procedures, they are in very good agree- 
ment and compatible within 2 <j error bars. Lastly, non- 



gaussianity tests based on higher order statistics were car- 
ried out, and do not show statistically significant departure 
from gaussianity at resolution 1 degree. The LGMCA and 
the official WMAP 9 maps essentially differ close to the 
galactic center where it remains extremely difficult to as- 
sess which map is less contaminated by foreground residuals 
or biases due to chance correlations in between CMB and 
foregrounds. 

The LGMCA code is available at 

http://www.cosmostat.org/lgmca and the LGMCA CMB 
map as well as the estimated power spectrum are available 
at http://www. cosmostat. org /product. 
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Appendix A: The LGMCA Method 

The GMCA Framework 

The GMCA (Generalized Morphological Component 
Analysis) method i s based on blind source separation (BSS) 
( |Bobin et al.| [2013). In the framework of BSS, each of the 
five WMAP frequency channels are modeled as a linear 
combination of n components: 



\fi = 1, 



i 5} %i 



(A.l) 



where Sj stands for the j-th component, &ij is a scalar that 
models for the contribution of the j-th component to chan- 
nel i and rii models the instrumental noise. This problem 
is more conveniently recast into the matrix formulation : 



X = AS + N 



(A.2) 



In practice, the number of components is set to n = 5 
which allows for more degrees of freedom to get a clean 
CMB map while keeping A invertible. Contrary to stan 



dard approaches in astrophysics (see (Bobin et al. 2013) 
and references therein), the GMCA relies on the sparsity 
of the components S in the wavelet domain. Taking the 
data to the wavelet representation only alters the statis- 
tical distribution of the data coefficients without affecting 
its information content. A wavelet transform tends to grab 
the informative coherence between pixels while averaging 
the noise contributions, thus enhancing the structure in the 
data. This allows to better distinguish components that do 
not share the same sparse distribution in the wavelet do- 
main. In addition, sparsity has the ability to be more sen- 
sitive to non-Gaussian processes, which has been shown to 
improve the foreground separation method. 

Having A as the mixing matrix and $ as a wavelet 
transform, we assume that each source Sj can be sparsely 
represented in 3>; Sj = aj<f>. The multichannel noiseless 
data Y can be written as 

Y = Aa$ . (A.3) 

where a is a N s x T matrix whose rows are ctj. 

This means the sparsity of the sources in 3? translates 
into sparsity of the multichannel data Y . The GMCA al- 
gorithm seeks an unmixing scheme through the estimation 
of A, which leads to the sparsest sources S. This is ex- 
pressed by the following optimization problem (written in 
the augmented Lagrangian form) 



1 



mm - 



IX- Aa$ 



Al 



(A.4) 



where typically p = (or its relaxed convex version with 

p = 1) and || X || F = sqrt (trace(X T X)^ is the Frobenius 
norm. 

Local GMCA 



The Local-GMCA (LGMCA) algorithm flBobin et al.|2013 ) 
has been introduced as an extension of GMCA: 

— multi-frequency instruments generally provide ob- 
servations that do not share the same resolution. 
For example, the WMAP frequency channels have a 



Band 


Obs. 


Patch size 


Res. 


I 


WMAP 9 data 
dust and Ha 


no 


60 


II 


WMAP 9 data 
dust 


256 


53 


III 


Ka,Q,V,W bands 
dust 


128 


39 


IV 


Q,V,W bands 
dust 


64 


30 


V 


V,W bands 
dust 


64 


22 


VI 


W band 
dust 


64 


15 



Table B.l. Parameters used in LGMCA to process the 
WMAP data with ancillary data. For each band, the second 
column gives the subset of data used to analyze the data, 
the third column provides the size of the square patches at 
the level of which the analysis is made, the fourth column 
gives the common resolution of the data. 



resolution that ranges from 13.2 arcmin for the W 
band to 52.8 arcmin for the K band. The makes the 
linear mixture model underlying the GMCA algorithm 
not valid. It is customary to alleviate this issue by 
degrading the frequency channels down to a common 
resolution prior to applying any component separation 
technique (the official CMB map provided by the 
WMAP consortium has a resolution of 1 degree). For 
this purpose, the data are decomposed in the wavelet 
domain and at each wavelet scales we only use the 
observations with invertible beams and then degrade 
the maps to a common resolution. This allows to 
estimate a CMB map with a resolution of 15 arcmin. 

most foreground emissions (e.g. thermal dust, syn- 
chrotron, free-free) have electromagnetic spectra that 
are not spatially constant. In the framework of GMCA, 
this translates into a mixing matrix A that also varies 
across pixel. Dealing with the variation across pixels of 
the electromagnetic spectrum of some of the compo- 
nents, the LGMCA estimates the mixing matrices on 
patches at various wavelet scales with band-dependent 
size. 



The LGMCA algorithm has been implemented and evalu- 
ated on simulated Planck data in ( Bobin et al.||2013 ). 



Appendix B: LGMCA parameters for WMAP data 

As described in Appendix [Aj LGMCA mixing matrices are 
estimated from a set of input channels at a given resolution, 
in a patch of data at a given wavelet scale. For WMAP data, 
the parameters used b y LG MCA t o com pute these matrices 
are described in Table |B.1| Figure [BT] displays the filters in 
spherical harmonics defining the wavelet bands where the 
derived weights (by inverting these mixing matrices) were 
applied to. 
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LGMCA filters in spherical harmonics 




Fig. B.l. Filters defining the wavelet bands used in 
LGMCA. 



Appendix C: Simulations 

In this section, the LGMCA algorithm is applied to data 
simulated by the Planck Sky Model ( PSM) developed by 
J. De labrouille and collaborator^] in ( Delabrouille et al. 



2012 ). The PSM models the astrophysical foregrounds in 
the range of frequencies probed by WMAP, the simulated 
instrumental noise and the beams. In detail, the simulations 
were obtained as follows. 

— Frequency channels: the simulated data are com- 
prised of the 5 WMAP channels at frequency 23, 33, 
41, 61 and 94 GHz. The frequency-dependent beams 
are perfectly isotropic PSF; their profiles have been ob- 
tained as the mean value of the beam transfer functions 
of at each frequency as provided by the WMAP consor- 
tium (9 years version). 

— Instrumental noise: instrumental noise has been gen- 
erated according to a Gaussian distribution with the 
covariance matrix provided by the WMAP consortium 
(9 years version). 

— Cosmic microwave background: the CMB map is 
a Gaussian random realization with theoretical power 
spectrum defined as the WMAP (9 years) best-fit power 
spectrum (from the 6 cosmological parameters model). 
The simulated CMB is perfectly Gaussian, and no non- 
Gaussianity (e.g. lensing, ISW, fjvz,) has been added. 
This will allow for non Gaussianity tests under the null 
assumption in the sequel. 

— Dust emissions: the galactic dust emissions is com- 
posed of two distinct dust emissions: thermal dust and 
spinning dust (a.k.a. anomalous microwave emission). 
Thermal dust is modeled with the Finkbeiner model 



((Finkbeiner et al. |1999[ )), which assumes that two 
hot/cold dust populations contribute to the signal in 
each pixel. The emission law of thermal dust varies 
across the sky. 

— Synchrotron emission: The synchrotron emission, 
as simulated by the PS M, is an extrapo lation of the 
Haslam 408MHz map ( flHaslam et al.|1982[ )). The emis- 
sion law of the synchrotron emission is an exact power 
law with a spatially varying spectral index. 

— Free-free emission: the spatial distribution of free- 
free emission is inspired by the Ha map built from the 



SHASSA and WHAM surveys. The emission law is a 
perfect power law with a fixed spectral index. 
— Point sources: infrared and radio sources were added 
based on existing catalogs at that time (including 
WMAP7 sources). In the following, the brightest point 
sources are masked prior to the evaluation results. 

A simulated WMAP dataset is produced for each of the 
9 years. This allows to process the simulated data in the 
exact similar manner as the WMAP data are processed. 

Component Separation 



2 For more details about the PSM, we invite the 



The same templates and parameters as listed in Table B.l 
were used for LGMCA. We also implemented an ILC as 
for the WMAP9 release: first computing the weights in the 
same regions as for the WMAP9 release, then smoothing 
them to 1.5 degree and finally applying them to the data 
at 1 degree in the same regions as defined in the official 
WMAP9 product. Note that no post-processing was per- 
formed to subtract the ILC bias due to foreground propa- 
gation as in the official product. This allows to compare on 
the simulations the relative performance of LGMCA and a 
localized ILC in pixel space at a resolution of 1 degree. 

Recovered Maps and Power Spectra 

The power-spectrum is comput ed fo llowing the procedure 
described in Section [3] Figure |C.1| displays the theoreti- 
cal power-spectrum in black and the LGMCA estimated 
power-spectrum in red. The pseudo-spectrum of the input 
map is shown in blue; these points would correspond to 
a perfect estimation of the CMB map where only cosmic 
variance is a source of uncertainty. The larger 1 o red er- 
rors originate from the error coming from the remaining 
instrumental noise. In this experiment, 75% of the sky cov- 
erage is used; the mask we used is a combination of point 
sources and galactic masks. These two plots show that the 
power-spectrum of the CMB map estimated after compo- 
nent separation does not show any statistically significant 
bias. 

The use of simulations allow for a precise decomposition of 
the CMB estimation error into its different components : 
CMB, remaining instrumental noise and foregrounds. For 
that purpose we apply the inversion parameters estimated 
with LGMCA independently to the simulated foregrounds 
and the instrumental noise. The resulting maps give the 
exact level of contamination of the CMB estimated by 
LGMCA. Figure [CT2] shows the power spectra of the CMB 
as well as the residual noise and foregrounds that contam- 
inate the estimated map. The resolution of the map esti- 
mated with LGMCA is 15 arcmin; therefore the different 
spectra in Figure |C.2| remains at the same resolution and 
are not deconvolved to infinite resolution. Again, exactly 
the same sky coverage of 75% is used in this experiment 
thus quantifying the exact level of foreground contamina- 
tion o f the estimated CMB power spec trum displayed in 
Figure [CTTj The two panels of Figure [C^2] first show that the 
main source of contamination is the remaining instrumental 
noise which predominates for i > 600. This translates to the 
large error bars of t he estimated power spectrum at small 
scales in Figure C.l For very low-^ (t < 20), the contribu- 



reader to visit the PSM website: \http://www. apc.univ- 
paris 7.fr/~ delabrou/PSM/psm. html 



tion of both the remaining noise and foregrounds is less than 
1% which is way below the error related to the cosmic vari- 



6 



J. Bobin et al.: WMAP 9-year CMB estimation using sparsity 



ance. In this experiment, the level of foreground contami- 
nation seem to be below 1% at all scales. This very low level 
has to be tempered : the ancill ary data, namely the com- 
posite all-sky H-alpha map of (Finkbeiner 200 31) and the 
Finkbeiner thermal dust template (Fin kbeiner et al7||1999 ) 
are also used within the PSM to produce the simulations 
of the free-free emission and thermal dust emission respec- 
tively. We should therefore expect the level of residual fore- 
grounds to be higher when LGMCA is applied to the real 
WMAP data. Interestingly, we also applied LGMCA with 
exactly the same parameters except that only the WMAP 
maps without ancillary data were used. The contamination 
levels are featured in Figure |C.2| in dashed line. If using 
templates indeed lowers the level of remaining foregrounds, 
their contribution is still much lower than the level of CMB. 
Noise is therefore the main source of contamination in the 
final CMB estimate whether external templates are used or 
not. 

Finally, in Figure |C.3| all components were propagated 
using weights computed by ILC and LGMCA. This figure 
illustrates that LGMCA is more efficient at lowering noise 
levels (due to the high amplification of noise in the ILC 
map when the 1 degree deconvolution is performed) and 
foregrounds contamination (due to localization and the use 
of templates) compared to the computed ILC map. 



Higher Order Statistics — Non-Gaussianities 

The level of non-Gaussianity in the recovered CMB map 
provides a sanity check to measure and localize any re- 
maining foreground contamination in the recovered CMB 
map, since in this case the CMB is generated as a Gaussian 
random field. In this work, we have computed the non- 
Gaussianity levels from the recovered CMB maps for 
LGMCA and ILC at 1 degree, and for the LGMCA map 
at 15 arcminutes. The 75 % mask was employed and sparse 
inpainting was performed to interpolate the signal inside 



the mask (Starck et al. 2013). The skewness and kurtosis 
were then computed on the simulation inside these masks 
on different wavelet sc ales using an isotropic un decimated 
wavelet on the sphere ( Sta rck, J.-L. et aT| 2006 ) , with the 
wavelet filters in spherical harmonic space described in 



Figure |C.4| These statistics were then centered on the ex- 
pected value (computed by propagating only the simulated 
noise and CMB) and normalized by the standard deviation 
computed from a set of 100 CMB and noise realizations. 
These statistics were also computed at different latitude 
bands for each wavelet scale to assess the level of foreground 
contamination in the maps at various scales and positions. 

The comparison between LGMCA and ILC at 1 degree 
is displayed in Figures |C.5| and |C.6[ For both methods, 
the skewness and kurtosis are compatible with the error 
bars due to propagated noise and cosmic variance, with 
a maximal detection at 2.5<r close to the galactic center. 
The same tests were also performed for the LGMCA map 
at the full resolution of 15 arcminutes and are displayed 



in Figures |C.7| C.8 and C.9 The difference observed be- 



tween the LGMCA non-gaussianity levels and those com- 
puted from the simulation without foregrounds is compat- 
ible with the errors expected at that resolution. 
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Fig. C.l. Estimated CMB map power spectrum from simulated WMAP (9 years). 
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Fig. C.2. Estimated CMB map, noise and remaining foregrounds power spectra from simulated WMAP (9 years) at 15 
arcmin resolution. 




Fig. C.3. Comparisons between LGMCA and ILC at 1 degree resolution. Estimated CMB map, noise and remaining 
foregrounds power spectra from simulated WMAP (9 years). Please note that the amplitudes have been further amplified 
by hj, the square of the 1 degree resolution Gaussian beam. 



8 



> 0.3 r 



J. Bobin et al.: WMAP 9-year CMB estimation using sparsity 
Wavelet filters for 60 arcminute maps Wavelet filters for 1 5 arcminute maps 




1000 



^ 0.2 



— Scale 1 

— Scale 2 
Scale 3 
Scale 4 
Scale 5 

— Scale 6 




Multipole 



Multipole 



Fig. C.4. Legendre coefficients of the wavelet filters employed for non-gaussianity analyses at 1 degree (left) and 15 
arcminutes (right). These wavelets are well localized in pixel space, allowing a fine analysis per latitude bands. 
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Fig. C.5. Comparison of skewness levels in LGMCA (red) and ILC (blue) maps at 1 degree computed for various wavelet 
scales. These statistics were centered on the expected value and normalized from a set of 100 simulations of CMB and 
Noise (see text). 
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Fig. C.6. Comparison of centered and normalized kurtosis in LGMCA (red) and ILC (blue) m aps a t 1 degree computed 
for various wavelet scales. The same mask and set of simulations where employed as in Figure [C3) 
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Fig. C.7. Centered and normalized skewness and kurtosis in LGMCA map at 15 arcminutes computed for various wavelet 
scales. A 75% mask and a set of 100 simulations of CMB and Noise were used to compute these statistics (see text). 



10 



J. Bobin et al.: WMAP 9-year CMB estimation using sparsity 



Wavelet Scale 1 Skewness 



Wavelet Scale 2 Skewness 



Wavelet Scale 3 Skewness 




Average Latitude (degree) 
Wavelet Scale 4 Skewness 





Average Latitude (degree) 
Wavelet Scale 5 Skewness 




:o o 50 

Average Latitude (degree) 

Wavelet Scale 6 Skewness 



Average Latitude (degree) 



Average Latitude (degree) 



Average Latitude (degree) 



Fig. C.8. Centered and normalized skewness in LGMCA map at 15 arcminutes computed for various wavelet scales and 
location. The same mask and set of simulations were employed to derive the statistic as in Figure [C]7 
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Fig. C.9. Centered and normalized kurtosis in LGMCA map at 15 arcminutes comput ed fo r various wavelet scales and 
location. The same mask and set of simulations were employed the statistic as in Figure [C?f| 



11 



